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1. INTRODUCTION 


Many of papers dealing wiih the spontaneous fission process 


describe it as a penetration of the nucleus through the potential 


barrier in a deformation space (1 A 2]k In order to have the pe- 
netrability (the barrier penetration probability ) P one usually 
iaatdes VKB formula 


B=[1+ exp (5445] ^ , (1) 
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where S43 is the reduced action given by the formule 


(2) 


$457 f N21v-815; da (2) 


(1) 


and it is taken along some trajectory given in prepared deforma- 
tion space. Here V = V(a) is the potential energy, Ba = B, (a) 
is the effective mass parameter (effective inertia) and E is 
the energy of the fissioning nucleus. The parameter "a" specifies 
the position of a point on the trajectory with a, and a, being 
the classical turning points corresponding to given energy E 
(V(a) = E). It might be the arc length of the trajectory consi- 
dered. 

The effective inertia EA associated with the fission motion 
along this trajectory is 


Mc Ma 
Bs B (a, — — G 
a > up 0 vocum ) 


where EM are components of the inertia tensor and dx 
(k21,2,...n) are deformation parameters specifying n-dimensional 
deformation space admitted in the description of the fission mo- 
tion. d 

The fission half life time Taf is calculated from the 
formula 


Ts = in 2 a (4) 


where n is the number of assaults of the nucleus on the fission 
barrier in time unit and is equal to the frequency eT of 
vibrations in the direction of fission degree of freedom. Further 
we shall assuming Rw, = 1 MeV inm zde 

There are two basic methods of selecting the path to fission 
in the case of multidimensional deformation space. 


i) The static method. The fission trajectory is taken to be 
that minimizing the potential energy Y((qj) [3]. 
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ii) The dynamic method. This method suggested originally in 
paper { 4) and developed for practical approaches in references 
f A Z 6] assumes that the fission path is that which mini- 
mizes the total action Sjo! This allows to include the influ- 
ence of Ba parameter on the process in a dynamic way. In ef- 
fect the fission path differs fron static one calculated in 
1) and the penetration probability 2 reaches its higher val- 
ue leading to the decrease of the spontaneous fission half life 
Tag [2 , 5). The minimization of S,5 might be realized in 
practice by applying Ritz method (see ref.[5 , 6]) or by dyna- 


mical programing methods (see ref.[2 , 5 , 6G). 


In both cases described above the penetrability $ is cal- 
culated along single path in deformation space. 

From the quantum mechanical point of view such à description 
of the fission process is incomplete. One may consider the fis- 
sion as a penetration of multidimensional fission barrier going 
through the all possible paths leading to a fragmentation of the 
nucleus. Each path should be taken properly with specific weight 
and the total penetrability should be a sum of all partial proba- 
bilities of fragmentation. Here the statistical mechanics comes 
into play and the method of averaging observables should be used. 
How to do all of that is the basic problem of our paper. 

In Section 2 the method of integration over the paths is 
proposed and the algorithm of calculating spontaneous fission 
half lives is described. 

Section 2 discuss the results of realistic model calculation. 
Different effects connected with the applied metropolis algorithm 
used in evaluating path integrals are considered, The results ob- 
tained in present work are compared with other ones derived on 
the basis of static 1) and dynamic ii) approaches (see eq. refs. 
[Lg deem ee em n 
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2. PATH INTEGRAL FORMULATION OF FISSION PROCESS 


2.1. Path integrals 


In this section we follow closely the way presented ín paper 
ff © 

Our basic expression is the Euclidean (imaginary time) ver- 
sion of Feynman path integral [7 , 8 , 9] 


z =fr e (5) 


The role of imaginary time plays an arc length "a" of the trajec- 
tory in n-dimensional (deformation) space {q} . The action s[q] 
is calculated on the actual path in deformation space and its 
form is given by eq. (2). 

Making the "a" axis discrete we adopt the following notation 


aaj) = 8 p ae OhnecOni (6) 


where - o¢ aX «9 , kzl,..,nD. 
The integral (5) reads then 


o y 
z= f n ag, exp(- $ S ((3]] > (7) 


The last expression is identical to the partition function 
for a statistical mechanics problem. The action S which couples 
nearest neighbour sites on a lattice Q;,_;» 3, and ARE plays 
the role of the hamiltonian of the considered system (e.g. a crys- 
tal). The Boltzmann factor is 


-1/5 S[3] 

e 
and the Planck constant h is equal to the temperature KT, 
The limit f —+» O picks out the classical configuration of the 
system(the single classical trajectory in our case). 
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Having the partition functicn it is simple to write average 
of any observable 0. In terms of the path integral one has 


rà -5/ f 
/ aq, O (qaos. eray )e 
(0) » 25.221 . ©; 
? PE | -s^ 
g dq, e 
-œ ísi 


Here dg, stays for the product dat »dq? s+. dare 


2.2. Averaged spontaneous fission half life time 


In the case of fission the quantity of interest 
taneous fission half life time Taf or shorter E. 
define as 


is ihe snon= 


f which ve 


(10) 


The quantity Ter appearing on the RHS of the formula (10) is 
the spontaneous fission hali life time as calculated in accord 
to the VKB prescription (4). The equation (10) is a statistical 
expression which contains all the informations about the nature 
of the potential energy V and the mass parameters B.  ( q ) 
entering the expression for the action S (eq. (2)). 3 
All global characteristics and the quantal features of the fis- 
Sion scenery are already contained in (1C). 

The average Tt contains contributions from ail possible 
trajectories to fission each of which contributes with proper 
weight determined by the action $[q]. 
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2.3. Monte Carlo tools 
The technique for evaluating the gum over trajectories is 
7 1 


based on the Monte Carlo method of Metropolis |7 
In this method an integral of the type 


[tū (1) 


is approximated by a sum o; a finite number of terns 


ce 


The points d in (10) are noi selecty! 2% fandom, but are dis- 
tributed in that region ot deżo 


rmatiiz Space saving the dominest 
contributions to the integral. In 


ase cl integrals given by 


Eqs. (8-10) the points 4, have io be selected ażcording tc the 
distribution 
exp ( - 8[41/% ] DQ 
Papa) = ———————— . (13) 


f Dą 2/8 


If such a selection of q, is done the Monte Carlo estimate O 
for the observable O simply reduces to an arithmetic average 
M 
1 
zg 2 Tee» (14) 
siet] 


o 
"m 


where M is the total number of configurations (trajectories) 
generated in the Monte Carlo sequence. 

The method of selection of important configurations is 
called "an importance sampling" [7 , 10 , 11]. Its realization 
is possible using a Markov process to generate M configura- 
tions (3 in equation (14). The process is constructed so 
that in the limit of large M the probability of occurence 
of a given configuration { à] is done by the Boltzmann distri- 
bution (15). 
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In the following ve sketch only the orastical rules of the 

polis algorithm [7] used in evaluating the Top Of eq. (10). 

Let EM jz1,...,1 Genote some initial (e.g. taken at 
"onüou) trajectory. Let us now select randomly a nen value for 
(for a specific  j) with uniform probability. Denote this new 


+4 


a. by "4.. If the action S is lowered by the replacement of 
a, t; 79, the variable at the site j is set to the new val- 
um "1 If however the change AS > O ‘then a random number 


r” with uniform probability distribution between O and 1 is 
generated and q, is changed to "a, only if 


exp (AS) yr’. (15 > 


In the last expression AS danotes the difference 
AS = sfa}- s(3] (16) 


where s[74] is the new action calculated along the new trajcc- 
tory passing through the point Ta and S is the action taken 
along the starting path. If the relation given by eq. (15) is not 
fulfilled the lattice variables at the site j retain their pre- 
vious values à E qj d). 

The next site (j+1) is probed in analogous manner and so on 
through ihe lattice, If this is made for all lattice points 
j21,2,...,N the one Monte Carlo step is done and a new trajectory 
i 3 ] is obtained, Repeating this process M times and then 
using eq. (10) the Monte Carlo average Th can be calculated. 

The Metropolis algorithm üescribed above is a one possible 
method of evaluating the path integral. Among other methods the 
"Heat Bath’ method is known 13) and "stochastic quantization’ 
method based on Langevin equation are in use especially in the 
field theoretical implementations of path integrals [15 , 16). 

Problems of convergence of the process of generating confi- 
gurations have been discussed widely in the literature [7 , 11 , 
12 , 15] and will be discussed in the concrete case of T.. cal- 
culation in the next section. 
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3. REALISTIC CALCULATIONS 


3.1. A model 

We consider a fission process in a 3-dimensional deformation 
stace -— A. ôn) where ZY is an average path to fission 
in the Aret (£5, £,) - space given by (2] 


E, = 0(0.05)1.00 
4 = 0.2 £, - 0.06 


The parameters às and 4, are pairing field parameters for 
protons and neutrons respectively. 

The potential energy V of the nucleus is obtained by the 
Strutinsky (18] prescription 


uidi xb TL. ue) 
where Erp is the liquid drop energy and the expression in 
brackets on the RHS is called a shall correction and it is a dif- 
ference between the BCS model energy and corresponding "'smooth^ 
energy of the nucleus (see e.g.; refs.[4 mmc ). In the fol- 
lowing we do not add a zero point energy to the potential energy 
surface [19] ). Mass parameters By, are obtained in the cranking 
approximation. The corresponding formulae is [3 bi 4] 


|<»|a8/ 3et|y'yf 


À (UyVy + Vy Wy) + Pige (19) 
Œy r Ryż 


By > d 2 
vv 


Here H is the single particle hamiltonian, uy and vy are 
the BCS occupancy parameters and Ey is the quasi particle energy 
corresponding to the single particle state |v). The term Pij 
describes the effect of the collective motion in the pairing 
field [4]. 

4 discussion of mass parameters coupled to the neutron and 
proton degrees Of freedom An and Ap is given in ref. [ 17]. 

The effective mass parameter as appearing in the action S 
(eq. (2)) reads in our case 
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- & E "a = ne 
= $ + 
Siege, a FE hi Aa" 
+ 2 £ A + B e A ; 
Će,, A, vV24 p " "t, Ân 724 p (20) 


where dot means the derivative taken over the path variable ‘a’. 
One can see that terms of the type B4 a are absent. 


n 
This is the effect of no mixing between —— and neutron pair- 
ing interactions which we assumed here. 


5.2. Monte Carlo experiment. Results and discussion 


In order to begin the Monte Carlo experiment we have to pre- 
pare the starting path to fission. From the calculations done 
early [? 3 525 11] it is known that curves of a physical signi- 
ficance are 
i) static trajectory (represented by a dashed line in Pig. 2), 
ii) dynamic trajectory (solid line in Fig. 2). 


LH - 
bale 


! log Te [yr] 


e Static start 


1 300 600 900 1200 1500 


Fig. 1. The convergence of the Metropolis algorithm. Full 

circles correspond to the logarithm of the averaged spon- 

taneous fission half life time 7T f as calculated after X 

Monte Carlo steps for the case of static starting path. 

Full triangles give the ln Tsf for the dynamic start.Error 
bars for both cases are also displayed. 
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This two different starts give us the possibility to check 
the convergence of the Metropolis algorithm. The convergence pro- 
cess is illustrated in Fig. 1. One can observe the change in the 
iogarithm of the mean value of Tar A5 calculated after ii Monte 
Carlo steps for both considered cases. On the beginning of the 
numerical experinent the values corresponding to static and dy- 
namic trajectories differ among themselves on about 6 units in 
log T... After M = 800 Monte Carlo steps (MCS) both logarithms 
of averaged Tof values coincide. The corresponding value of 
log T,, is equal to 4.70. This result changes to the value 4.40 
in the successive 700 MCS. The overlap reached gives an evidence 
thet the Monte Carlo process applied is fairly convergent. Hov- 
ever the increasing of EX may lead to the unexpected jumps in 
ine log Tor values which increase slightly the result 4.40. Such 
a behaviour cf the Monte Carlo process is a consequence of the 
appearance of some "wila" trajectories on which the action S is 

ery large one. This leads in consequence to the increase of 
average Tore The influence of such undesired curves can be eli- 
minated after a very long runs of the lente Carlo program. On the 
otner hand in long runs the periodic nature of the pseudo random 


: 


Na east, anm t 


Fig. 2. Dynamic (solid line) and statica (dashed line) paths 
to fission an WEDS TA y Ay) ‘space for ?2?7pg, Pull circles 
rcpresent the avêfage” path to fission as calculated after 
1500 Monte Carlo steps starting from the static path. Full 
triangles represent the average trajectory for which the 
starting curve was the dynamic path to fission. 
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number generators may cause the additional complications. There- 
fore ve have stopped the Honte Carlo process at 1500 steps. 

All results presented here are obtained for 252m nucleus 
for which the experimentally known value of log Taf is equal to 
2 and our predicted value differs on two units from it. At the 
same time the Monte Carlo result is smaller than the static value 
of log Toe’ 2 5.5. 

Fig. 2 shows starting paths to fission: static (dashed line; 
and dynamic (solid line). Full circles and full triangles show the 
average trajectories obtained after 3! = 1500 MCS starting from 
static and dynamic curves respectively. Averaged trajectories in 
( fog: 5 y! An) space differ insignificantly. The corresponding 
standard deviations are also shorm. 


> ZZ LL. 


0.2 04 0.6 0.8 
Fig. 5. Dynamic (Coe toe] and static (dashed line) 


fission barriers for Fm as compared to averaged bar- 

riers calculated after 1500 Monte Carlo steps. Full cir- 

cles and triangles correspond to static and dynamic 

starting trajectories respectively. Statistical error 
bars are also shown. 


In Fig. 3 corresponding fission barriers are displayed as a 
function of elongation coordinate E One can see that both 
averaged barriers agree well within an accuracy given by the stan- 


dard deviations. Averaged tarriers are situated between static and 
dynamic barriers. 


A. Baran, A. Staszczak 


Fig. 4. Effective mass parameters Bip, as calculated 


along static (open circles), dynamic (open triangles) and 
Monte Carlo averaged paths: static (full circles) and dy- 
namic (full triangles). 


The next Fig. 4 shows the effective mass parameter Barr 
calculated along the average paths to fission (given in Fig. 2) 
as compared to masses taken along static and dynamic paths to 
fission, Again, full circles correspond to the static start of 
Monte Carlo process and full triangles correspond to the dynamic 
start. In opposition to barriers, effective masses were calcula- 
ted on average paths and they are not averaged Monte Carlo messes, 

It is too early to conclude that the present path integral 
approach to spontaneous fission process gives results which are 
better then that reported early [2-5, 17]. However, statistical 
model of the process may lead to the more correct insight into 
the fission fenomenon and may indicate the way for many dimen- 
sional analysis of it. 

The investigations cornected to a systematics of spontaneous 
fission half life times, hights of fission barriers, an asymptotic 
behaviours of effective masses etc. are in progress. 


Authors are very thankfull to all of their colleagues who 
were and who are patient of discussions about the Monte Carlo 
stuff and path integrals. 
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STRESZCZENIE 


W pracy zaproponowano metodę obliczenia połówkozych czasów 
życia jąder atomowych ze względu na rozszczepienie spontani czne 
stosując formalizm całek po drogach i mechanikę statystyczną, 
Obliczenia przeprowadzono w trójwymiarowej przestrzeni parametrów 
Op An oraz 6,» gdzie A jest przerwą energetyczną typu BCS, 
a >; jest parametrem deformacji przestrzennej jądra, 

Otrzymane czasy życia są większe njż te, które obliczano 
wcześniej stosując minimalizację całkowitego działania i biorąc 
w rezultacie tylko jedną drogę do rozszczepienia, 


PESHME 


B paóore npeiuoraeTCs METOM BHUKCJEHAA BpeMeH RUSHM ATOMHNX 
HZep MO OTHOZeHMD K CMOHTAHHOMY 1EN6HKO C NpHMEHGHMGM QOpMa- 
USHA MHTerpasos NO TPAEKTOPMAM KH CTATMCTUYECKOŃ MexaHHkKH. Pac- 
YeTH npoBOAMIZMCP B TpexMepHOM IpocTpaHCTBe HapaHeTpOB Ap, An 
u fous rae A - 3T0 2Heprerzyeckag wenb Tuna BKH, a bon - na- 
paueTp npocTpaHcTseHHOÁR AZedopuanAut AAPA. 

TlonyyeHHNe npeueHa KUSHK ÓONBWe Tex 3HaueHul, KOTODHE BH- 
WHCARAMCE paHbue C NOMOLBO MKHKMK3ANKK Nonsoro AeKcTBMA M pac- 


CMOTpeHMH TOHBKO OĄHOTO KaHaJa ACACHUA. 


